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Abstract 

The most general one dimensional reaction-diffusion model with nearest- 
neighbor interactions, which is exactly-solvable through the empty interval 
method, has been introduced. Assuming translationally-invariant initial 
conditions, the probability that n consecutive sites are empty (E n ), has 
been exactly obtained. In the thermodynamic limit, the large-time behav- 
ior of the system has also been investigated. Releasing the translational 
invariance of the initial conditions, the evolution equation for the proba- 
bility that n consecutive sites, starting from the site k, are empty (Ek,n) 
is obtained. In the thermodynamic limit, the large time behavior of the 
system is also considered. Finally, the continuum limit of the model is 
considered, and the empty-interval probability function is obtained. 
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1 Introduction 



The principles of equilibrium statistical mechanics are well established. But 
thermal equilibrium is a special case, and little is known about the properties 
of systems not in equilibrium, for example about the relaxation toward the sta- 
tionary state. There is no general approach to systems far from equilibrium. As 
mean-field techniques, generally, do not give correct results for low-dimensional 
systems, people are motivated to study exactly-solvable stochastic models in 
low dimensions. Moreover, solving one-dimensional systems should in principle 
be easier. Different methods have been used to study these models, including 
analytical and asymptotic methods, mean-field methods, and large-scale numer- 
ical methods. Exact results for some models on a one-dimensional lattice have 
been obtained, for example in [Q-^6|. 

The term exactly-solvable have been used with different meanings. In fT^| , 
a ten-parameter family of reaction-diffusion processes was introduced for which 
the evolution equation of n-point functions contains only n- or less- point func- 
tions. The average particle-number in each site has been obtained exactly for 
these models. In Jl9| , the same method has been used to analyze the above men- 
tioned ten-parameter family model on a finite lattice with boundaries. In p8| 
and Jig] , integrebility means that the iV-particle conditional probabilities' S- 
matrix is factorized into a product of 2-particle S-matrices. 

The empty interval method (EIM) has been used to analyze the one dimen- 
sional dynamics of diffusion-limited coalescence p0|-p3[. Using this method, 
the probability that n consecutive sites are empty has been calculated. This 
method has been used to study a reaction-diffusion process with three-site in- 
teractions [^4| . EIM has been also generalized to study the kinetics of the g-state 
one-dimensional Potts model in the zero-temperature limit |2a ] . 

In this article, we are going to study all the one dimensional reaction-diffusion 
models with nearest neighbor interactions which can be exactly solved by EIM. 
It is worth noting that ben-Avraham et al. have been studied one-dimensional 
diffusion-limited processes through EIM |2p|-^]. In their study, some of the 
reaction rates have been taken infinite, and they have worked out the models on 
continuum. For the cases of finite reaction-rates, some approximate solutions 
have been obtained. 

We study models with finite reaction rates, obtain conditions for the system 
to be solvable via EIM, and then solve the equations of EIM. We do this for a 
system on a lattice and on continuum. 

The scheme of the paper is as follows. In section 2, the most general one 
dimensional reaction-diffusion model with nearest-neighbor interactions which 
can be solved exactly through EIM has been introduced. Assuming translational 
invariance, the probability that n consecutive sites are empty (E n ), has been 
exactly obtained. In the thermodynamic limit, the large-time behavior of the 
system has also been investigated. In section 3, the assumption of translational 
invariance has been released, and the evolution equation for the probability that 
n consecutive sites, starting from the site k, are empty (Ek,n) is obtained. In the 
thermodynamic limit, the large time behavior of the system is also considered. It 
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is shown that translationally-asymmctric fluctuations relative to the stationary 
configuration disappear faster than the translationally-symmetric fluctuations. 
In section 4, the continuum limit of the model is considered, and the empty- 
interval probability function is obtained. 

2 Models solvable through the empty interval 
method: the translationally-invariant case 

Consider a general one-species reaction-diffusion model on a one-dimensional 
lattice with L + l sites, with nearest-neighbor interactions. We want to impose 
restrictions on the reaction- and diffusion-rates, so that the system is solvable 
via EIM, that is, so that the evolution equation for the probability that n 
consecutive sites are empty (E n ) is closed. Suppose that the initial condition of 
the system is translationally-invariant. The most general interactions for a single 
species model in a one-dimensional lattice with nearest-neighbor interactions are 

• o — > (••, o», oo) o • — ► (»o, ••, oo) 

•• — ► (»o, oo, o») o o — > (»o, ••, o»), (1) 

where an empty (occupied) site is denoted by o (•). The constraint of solvability 
of the model through EIM, imposes that, as we will show, there are no processes 
in which the final configuration is oo. We shall also see that the processes the 
initial configuration of them is oo, have no effect on the solvability through EIM. 
But first, let us consider only the systems for them there are no interactions with 
oo as the initial or final configuration. So, among the above 12 interactions, only 
the following 6 interactions remain to be considered. 



I ••, n I »o, r 3 

The parameters r, are the rate of interactions. Define 




(2) 



PfT^S) =: E n , (3) 

from which, one obtains 

n n 

P(m ZZ^TS) = P(Z^S •) = E n -E n+1 , (4) 

where • (o) indicates an occupied (empty) site, and P denotes the probability 
of the configuration. The evolution equation for E n (t) is 

dE n n 

=r 5 P(»To*~~o) + r 3 P(o9~o^~o) 

n n 

+ r e P(o • • • o • •) + r 2 P(o • • • o • o) 

n n 

- (ri + r 2 )P{» o~o*~~o) - (r 3 + r 4 )P(o^*o~o •). (5) 
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The right-hand side of (|J) is not generally in terms of only E n 's. In order to 
make it expressible in terms of E^s, using (Q), one has to impose 

^3 = r 5 , r 2 = r 6 . (6) 

Then, using 

n n n 

P(.7^S) + P(o7^~$) = P(7^S) = E n -! - E n , (7) 
and another similar relation, one arrives at 

= (ra + r 3 )(£„_i + E n+1 - 2E n ) - (n + n){E n - E n+ i), n > 1. 

(8) 

Note that if there were reactions with the final configuration oo then one en- 
counters with terms like P(o o • o • • • o) or P(o o • • o • • • o) at the right-hand 
side of (|5|), which are not expressible in terms of Ek's. On the other hand, if oo 
is the initial configuration in a reaction, the evolution equation for E n 's is still 
closed, although the analogue of (|J) will be a linear finite difference equation 
with nonconstant coefficients. 

The equation of motion of E\ (t) is 

= (r 2 + r 8 )(l + E 2 - 2E 1 ) - (n + r 4 )(^ - E 2 ). (9) 

at 

It is seen that it takes a form similar to (|§|) , provided one defines 

E (t) := 1. (10) 

Then we have (||), for n > 1, equipped with the boundary condition ([To]). We 
also set £x + i(t) = 0, which means that initially at least one particle is present 
in the lattice. If initially all the sites were empty (E n (0) = 1 for all n), then the 
above-defined reactions would not change the configuration of the system, and 
if initially at least one particle was present, then at any time EL + \(t) = 0. So, 
the completely empty lattice is a stationary state which is decoupled from any 
other state. Defining 



(11) 



7 n + r 4 
b ■= — - — , 

r 2 + r 3 

and rescaling the time properly, the equation of motion becomes 

=E n _ 1 + E n+1 - 2E n + b(E n+1 - E n ), 0<n<L+l (12) 

at 

with the boundary conditions 

E (t) = l, E L+1 (t) = 0. (13) 
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A particular solution to this is the stationary solution: 

<-i + El +1 - 2El + b{El +l - El) = 0. (14) 

Taking the ansatz 

E%=Az?+Bz$, (15) 

for Ef i: one arrives at 

z l + z7 1 -2 + 6(^-1) =0, (16) 

the solutions of which are z\ — 1/(1 + 6) and Z2 = 1. Using the boundary 
conditions Eq = 1 and -Ex+i = 0, A and i? are obtained as 

A 



1- (1 + 6) 



-i-i 



-(l + fe)- L - J 
l-(l + 6)" 



^ r^^U-r (17) 



Defining 

F„(t) := £„(*) - El, (18) 

it is seen that the evolution equation for F n is the same as that of E n , but the 
boundary conditions for F n are homogeneous. The initial condition for F n is 

F n (0) = E n (0) - El. (19) 

To calculate F n (t), one seeks the eigenvalues and eigenvectors of the operator 
at the right-hand side of (|l^), that is 

e fn = fn+l + fn-1 ~ Vn + &[/n+l - fn}- (20) 

The solution to this is 

f n =az? + bz$, (21) 

where z,'s satisfy 

z J 2 (l + 6)-z J (e + 2 + 6) + l = 0. (22) 

Now, defining 



:= z,Vl + 6, (23) 

it is seen Z1Z2 — 1. So 

/» = 77^77- M n + in (24) 

(1 + 0) 2 
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The boundary conditions Fo(t) — i<L+i(i) = 0, lead to 



where k is an integer satisfying 1 < k < L + 1 , and 



Z = CX P 777 > ( 25 ) 



/fe '" = (TT^ sin (^ ) " ,2< " 



The corresponding eigenvalue is then 



e k = -2-b + 2VTTfocos ( ) . (27) 



A: 



otk . I imk 



Then F n (t) will be 

^g^^K*' (28) 

where 

^ = ITT E ^(°) - + ^)™ /2 (^) • (29) 

m— 1 ^ ' 

In the thermodynamic limit (i — > oo), F n (t) takes a simpler form. Defining 
x := nk/(L + 1), (|f) and (||) lead to 

^n(«) =- V r dx e ( - 2 - b+2vmcosa; ) t sm(na ; )sin(TOa ; )F m (0)(l + fe) (m - ,l)/2 

OO 

= (1 + &) (m -" )/2 e- (2+&)t (I m -„(2iVTT5) - l m+n (2tVTTb)) 

m=l 

x[E n (0)-E*], (30) 

where in the second line we have used the integral representation of the modified 
Bessel functions 

!*(*) = -/" da;cos(nx)e tcos;r . (31) 
* Jo 

To study the large-time behaviour of F n (t), one takes A k =: nBk/(L + 1). 
In the thermodynamic limit (£ — * oo), ( p8| ) leads to 

= (TtW if d ^W e(_2 ^ +2yTT5cOSa;)t sin (^)- (32) 
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At large times, the main contribution to the integral comes from the region 
x rs 0, in which the exponent of the exponential term takes its largest value. 
So, 



or 

(-2-b+2 v / T+fc)t 

F " (<) i ' (34) 

provided B(x) is well-behaved and nonzero at x = 0. If 6 ^ 0, there is an 
energy gap in the spectrum and the system relaxes towards its stationary state 
exponentially If b = 0, there is no energy gap and the relaxation towards the 
stationary state is in the form of power law with the exponent — 1. 

The empty-interval probability functions can be used to obtain some kinds 
of n-point functions. It is easy to see that 

tn 

P(.^T75.) = E m ~ 2E m+1 + E m+2 . (35) 

So, the results for E n (t)'s can be used to obtain the probability that between 
two occupied sites, there are n sites, which are empty. 



3 Models solvable through the empty interval 
method: the general case 

In the previous section, we considered translationally-invariant initial condi- 
tions. As the dynamics is translationally-invariant, the probability E n {t) will 
be the same for all sites, provided the initial condition for it is so. In this section, 
we release the translational invariance of the initial conditions and the quantity 
of our interest is the probability E ktTl (t), that n consecutive sites, starting from 
the site k are empty at time t: 



P k (^-~5) = E k , n (36) 



It is easy to see that 



F;(»oo • • -o) = Ei + i jm - E, 



P;(oo ■ • -om) = Ei, m - Ei >m+1 - (37) 
Using the interactions (^|), and the above identities, one arrives at 

^~ ' ='"3(-E'fe+l,n-l — Ek,n) + r 2(Ek,n-l ~ 

- (n + r 2 )(E kiTl - E k ^x,n+x) - (r3 + ri){E k , n - E kn+1 ). (38) 
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Similar to the previous section, imposing the boundary condition 

E kfi (t) = 1, (39) 

makes the evolution equation valid for < n < L + 1. If the lattice is not 
initially empty one also has 



E k , L+1 {t) = 0. (40) 



Using the definitions 







r-2 


+ r 3 




+ r 2 


r-2 


+ r 3 




r 3 


'2 


+ r 3 



d:=^—, (41) 
equation (|3^) can be rearranged in the form 

k ^ - =-Efc,n-l + Ek : n+1 — 2Ek, n — b(Ek, n — Ek,n+l) 

— c{Ek, n +i — Ek-\,n+l) + d(Ek+l^ n -l — Ek, n -l)- (42) 

Using the particular solution E%, one defines 

F k , n (t) := E k , n (t) - E*, (43) 
which satisfies (ff2|), but with homogeneous boundary conditions: 

F k ,o{t) = F krL+1 (t) = 0. (44) 
Applying the Fourier transformation 

F n (w,t) :=5^w*F fc)Tl (t), (45) 



one arrives at 
dF n (w,i) 



F w _i(w, t) + F„ +1 (w, t) - 2F„ (w, t) - 6[F„ («, t) - F n+1 (w, i)] 



di 

- c[F n+ i(w,t) - w# n+ i(w,t)] + d[a; _1 jF„_i(w,t) - F n _i(w,i)]. 

(46) 

To solve this, one first solves the eigenvalue problem 

e/nM =/ n _i(u>) + /n+lM - 2f n (uj) - b[f n ((j) - f n+1 (uj)} 

~ c[/ n +l(w) - w/„+i(w)] + d[w _1 /n-lM - /n-lM]. (47) 
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using the ansatz 

f n = az? + bz2, (48) 
and the boundary conditions, one arrives at 



f s , n (u) = [B{w)] n l 2 sin 



nirs 



L + l / 

2-& + £)(w)cos( T ^f T ), (49) 



where 



D(u>) :=2 V[l + diuj- 1 - 1)][1 + b + c{lu - 1)], 

v ; l + 6 + c(w-l) v ; 

and s is an integer between 1 and L. In the thermodynamic limit L — > oo, one 
arrives at 

F n (w, i) = £ F m (w, 0)[i?(c)]("- m )/ 2 e -( 2 + b )* 

m 

X{I 

rn — n [D(oj)t]-l m +n[D(oj)t]}, (51) 



where 



F m (cj,0)=J2^ k [Ek, m (0)-E^}. (52) 



Now, let's consider the relaxation of the system towards its stationary state. 
Suppose that the initial value for Ek tU is so that F n {oj, 0) contains a term pro- 
portional to S(p) (where u) = e w ), and another term which is a smooth function 
of lo. The delta term comes from a translationally-invariant part in F n ^(0). 
Using the steepest descent method, one can see that the relaxation behavior of 
the second term is governed by the extremum value of the eigenvalues e with 
respect to a complex lo. This is found to be 

£max = -2 - b + 2[Vd~c +y /(i-d)(l + b- c)]. (53) 

It is easy to show 



Equality holds when 



< -2-6 + 2V1 + 6. (54) 



r2{n +r 2 ) = r 3 (r 3 +r 4 ). (55) 



This means that the relaxation time for the translationally-noninvariant part is 
smaller than of the translationally-invariant part. That is, the translationally- 
noninvariant fluctuations disappear faster than the translationally invariant 
parts. 
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4 The continuum limit 

In the previous section, we considered the probability of finding n consecutive 
empty sites starting from the fc-th site. For the continuum limit, it is better to 
use a quantity with arguments symmetric relative to the starting point and the 
end point of the empty interval. For an empty interval of the length n starting 
from the site k, the end site is k' = 2k + n — 1. Then one can use 

s := k + k' = 2k + n- 1 (56) 

instead of k for labeling the empty interval, s/2 is the center of the empty 
interval. So, one uses the quantity 

SsAt) ■= E ktn (t). (57) 
The equation of motion for T s ^ n {t) is then 

^ - - ='"2(^ r s-l,ri-l +-7 r s-l,n+l _ 2jF s rl ) + r$ (Fs+l ,n-l + ^s+l.n+l — 2T St n) 

+ ri(T s+ i, n+ i - T s , n ) + ri(T s -i. n+ i - T s , n ). (58) 

Here T is the solution to the evolution equation of £, but with homogeneous 
boundary conditions. Using X := s/2 and x :— n, and Taylor-expanding the 
above expression in the continuum limit, one arrives at 

d^ 3 ^ = (Ad x + Bd x + ^d 2 x + Cd 2 x + Dd x d x )T(X, x; t), (59) 
where the parameters A, B, C, and D are 

A := r 3 - r 2 + 4 1 , B := n + r 4 
C:=±[r 1+ r 4 + 2(r 2 +r 3 )], D:= T -±^. (60) 
Using the change of variables 

±:=X + {A -^ )t -2^ X > (61) 



and 



one arrives at 



B B 2 

T(X, x- 1) =: cxp[-— x - -git]T{X, x; t), (62) 



dHX dt X]t) - \c% + ( j - *; *)■ (63) 

The boundary conditions for T are 

fi(X,x = 0;t)=fi(X,x^oo;t)=0. (64) 
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The Green function G(X , X' , x, x';t) for the equation ( |63| ) with the boundary 
conditions (pi) is 



G(X, X', x, x'; t) := l e -(*-*'V(4C't) 



X 



e -(x-x')/(4Ct) _ e -(x+x')/(4Ci) 



(65) 



where C" := (C/4) - L> 2 /(4C). Finally, 



:F(X..r:l) = exp j 

OO fOO 



pOO poo 

/ / Ax'dX' GiX.X^x.x'-^^TiX'.x';^ 

Jx'=Q J X' = -oo 



/x'=0 J X' = -oa 

(66) 



To obtain the solution for £, one has to add a particular solution with the 
boundary conditions 

£(x = 0) = l, £{x^oo) = 0. (67) 

One can choose this particular solution to be translationally-invariant (that is 
X-independent). One is then led to 

(Bd x +Cag)£ p (x)=0, (68) 

the solution to which is 

£ p = cxp(-Bx/C) (69) 
This particular solution is the same as that of section 2 in the continuum limit. 
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